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I ' Using the density matrix renormalization group technique, we study the ground 

^ ■ 

o 
o 



state phase diagram and other low-energy properties of an isotropic antiferromagnetic 
spin-half chain with both dimerization and frustration, i.e., an alternation S of the 



^ ■ nearest neighbor exchanges and a next-nearest-neighbor exchange J2 . For (5 = 0, the 



system is gapless for J2 < J2c and has a gap for J2 > J2c where J2c is about 0.241. 
For J2 = J2c, the gap above the ground state grows as S to the power 0.667 ±0.001. In 
the J2 — S plane, there is a disorder line 2J2 + 6=1. To the left of this line, the peak 
in the static structure factor S{q) is at Qmax = (Neel phase), while to the right of the 
line, qmax decreases from tt to 7r/2 as J2 is increased to large values (spiral phase). For 
(5=1, the system is equivalent to two coupled chains as on a ladder and it is gapped 
for all values of the interchain coupling. 

PACS numbers: 75.10.Jm, 75.50.Ee 
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I. INTRODUCTION 



One of the most widely studied spin models is the isotropic spin-half Heisenberg 
antiferromagnetic chain. This model is Bethe ansatz soluble^'^ for uniform nearest- 
neighbor (NN) exchange constants. Since we can map the spin-half model in one 
dimension to a system of interacting spinless fermions by the Jordan- Wigner transfor- 
mation, we expect a Peierls type of dimerization in this system when it is coupled to 
phonons. The dimerized model, with exchange constants alternating as 1 ± 5, is not 
exactly solvable, although many numerical and approximate analytic results exist for 
this model. There has also been considerable interest in the spin-half chain^ with both 
NN and next-nearest- neighbor (NNN) antiferromagnetic interactions J2. Interestingly, 
the ground state of this model can be solved exactly^ for J2 = 0.5. The ground state is 
doubly degenerate for the ring and the two states are representable as the two Kekule 
structures of the valence-bond theory. 

The heuristic phase diagram (Fig. 1) of the ground state for a model that incor- 
porates the effects of both dimerization and frustration, i.e., the J2 — S model is known, 
although there has been no systematic study of this model in the J2 — S plane. The 
ground state is exactly known for points on the line 2J2 + 6 — 1 ^. Analytical and nu- 
merical studies^"^ of the model on the S—0 line show a transition from a gapless phase 
for J2 < J2C to a gapped phase for J2 > J2C ■ The value of J2c has been accurately 
computed to be 0.2411 ± 0.0001 in Ref. 8. The J2 = 0,6 > line corresponds to a 
dimerized spin chain. Another interesting line on the phase diagram is 5 = 1 which 
corresponds (cf. Fig. 2) to coupled spin chains with J2 as the intrachain and 2 as 
the interchain coupling constants. Little else is known about the general J2 — S model. 
Numerical studies of this model have been largely confined to exact low-lying states of 
small systems A?" < 22 on the 5 — line. 

In this paper we present a systematic study of the various regions of this phase 
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diagram using the density matrix renormalization group (DMRG) method^. Using this 
technique we have studied both open chains and rings of sizes ranging between 200 and 
300 sites, depending upon the parameters of the model. In section II, we give a brief 
introduction to the method and present some computational details. In section III, we 
present our results and discuss them. 

II. THE HAMILTONIAN AND THE DMRG TECHNIQUE 

The Hamiltonian for the spin-half chain we have studied is given by 

H = Hq + Hi + H2 

2N-1 

Hq — ^ Si ■ Si^i 
i=l 

2N-1 

Hi ^ - s J2 i-'^ySi-Si+i 

i=l 
2N-2 



H2 — J2 ' 1 



i=l 



where J2 is the NNN exchange and S is the NN exchange alternation parameter. (Note 
that we have set the average NN exchange Ji — 1). In the open chain with an even 
number of sites, the dimerization is chosen to make the exchange constant for the 
terminal bonds strong (1 + 6). The topology of the chain being studied is shown in 
Fig. 2. 

The DMRG technique involves systematically building up the chain to a desired 
number of sites starting from a very short chain by adding two sites at a time. The 
initial chain of 2n sites, with n a small enough integer, is diagonalized exactly and 
the reduced density matrix for the left n sites is computed from the ground state of 
the 2n chain Hamiltonian by integrating over the states of the right n sites. The 
density matrix is diagonalized and a matrix representation of the n-site Hamiltonian is 
obtained in a truncated basis using m basis vectors which are the eigenvectors of the 



density matrix corresponding to its m largest eigenvalues. The Hamiltonian matrix 
for the 2n + 2 chain is then obtained in the (2s + l)^m^ dimensional direct product 
subspace obtained using the truncated basis of the left and the right parts of the 2n 
chain and the full space of the two additional spins which are inserted in the middle. 
After obtaining the ground state of the 2n + 2 chain in the truncated basis, the density 
matrix of half the chain, now with n + 1 sites, is obtained. The procedure is repeated 
to obtain the ground state of the 2n + 4 chain. The iteration is stopped on reaching 
the desired chain length. 

The DMRG technique works best for spin-half chains with short-range interactions. 
The accuracy depends crucially on the number of eigenvalues of the density matrix, 
m, that are retained. For spin-half chains with NN interactions only, m = 32 gives 
results that approximate very well to the infinite system results. In systems with 
longer range interactions, it becomes necessary to keep track of the spin matrices in 
the most recent basis, corresponding to all the sites which interact with the new spins 
that are introduced at a given DMRG step. While the DMRG algorithm is devised 
to minimize the errors due to truncation of the Hamiltonian matrix, the algorithm 
does not necessarily retain the spin matrices very accurately. The limited success 
of the DMRG technique for two-dimensional lattices can be attributed to this fact, 
since a two-dimensional lattice is topologically equivalent to a one-dimensional lattice 
with very long-range interactions. In the case of spin-half chains with NN and NNN 
interactions, the spin matrices that appear in the Hamiltonian at any stage would have 
undergone at most two transformations, and the results we find for m=64 compare well 
with exact numerical diagonalizations of chains upto 22 sites. 

III. A FIELD THEORY FOR THE NEEL PHASE 

Before presenting our numerical results, we briefly discuss a fleld theory for the 
low-energy and long-wavelength modes of the spin chain. In the Neel phase, the low- 
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energy physics is well described by an 0(3) nonlinear sigma model (NLSM) with a 
topological term^'^'^^. (The field variable is an unit vector). This led Haldane to predict 
that integer and half-integer spin chains will have qualitatively different behaviour for 
5 = 0, in that the former have a gap in their excitation spectrum while the latter do 
not. This prediction has received support from numerical calculations and experiments. 

Since our primary interest in this paper lies in the spin-half chain, it is convenient 
to use a different approach which is specific to these chains. This is the technique of 
bosonization discussed by Luther and Peschel^^ and AfHeck^^. In this approach, we 
first carry out a Jordan- Wigner transformation on the spin chain to obtain a model of 
interacting spinless fermions on a lattice 

St = V'i exp [ iTT ^ ipjipj ] . ^ 

j<i 

The spin-up and spin-down states at a site correspond respectively to states with 
and without a fermion. The continuum limit of the transformed Hamiltonian can be 
bosonized to obtain a tractable bosonic theory in 1 + 1 dimensions. For our system, 
this just produces the theory of a free 'relativistic' massless boson. The corresponding 
Lagrangian density is given by 

c ^ I {d,<t>f , (3) 

where is a bosonic field. The scaling dimension of an operator exp(iQ!0) is therefore 
c? I^.'K. The relation between a spin operator and the continuum field is given by 

S\x) - sin(V2^0) , (4) 

where a is the underlying lattice constant. Since the spin Hamiltonian is isotropic, 
it is sufficient to consider the S^S^ correlations. A straightforward evaluation of this 
correlation function yields 

( S\x) S%0) ) ^ ^--l— . (5) 
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This power-law correlation implies the absence of a gap above the ground state. Eq. (5) 
generally has corrections like log(a;/a) because the field theory has a marginal operator 
O. In the spin chain, this corresponds to the NNN exchange J2 . 

There is a crtical value J2c such that O is marginally irrelevant for J2 < J2c 
and marginally relevant for J2 > J2C • At J2C , the field theory is exactly conformally 
invariant and there are no logarithmic corrections. The spectrum at that point is 
described by a specific conformal field theory. In particular, the gap between the 
ground state and the low-lying excited states scales with the chain length as 1/N. 
Further, the first excited states with S = and S = 1 are exactly degenerate. For 
J2 < J2c 5 the excited S — ^ state has a higher energy than the S — 1 states. For 
J2 > J^c 5 the 5" = state becomes degenerate with the ground state as A?" — > 00, and 
the triplet states are separated from the two ground states by a finite gap. One can 
therefore determine J2C accurately as the point where the singlet and triplet excited 
states are degenerate for large values of N'^'^. 

We can treat the spin-Peierls term Hi in Eq. (1) as a perturbation on Hq for 
small values of 5. Using the above procedure, we obtain the bosonic representation of 
Hi in the continuum limit as 

Hi ^ 5 cos(\/2^ (j)) . (6) 
This is a relevant operator with scaling dimension i and it therefore produces a gap 

A ^ 52/3 . (7) 

Similarly, the change in the ground state energy caused by the perturbation in (6) 

scales as 

Eo{5) - £^o(0) - 5^/3 . (8) 

These results are perturbative^^ and they are valid only for small 5. Further, there are 
generally corrections of order log ( 5 ) in Eqs. (7) and (8) because of the presence of 
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the marginal operator. However if we are exactly at the critical point J2c , there are 
no logarithmic corrections to (7) and (8). 

For J2 > J2C and 5 = 0, there is a gap between the two degenerate ground states 
(both S — 0) and the first excited state (S — 1). Near , the gap A has an essential 
singularity of the form 

A - exp ( - — ) (9) 

Hence the gap is numerically indistinguishable from zero unless J2 is greater than 
about 0.3. 

Before ending this section, we note that a different field theory is required in the 
spiral phase (which, for (5 = 0, sets in beyond J2 = 0.5 for spin-half). This field theory is 
a NLSM based on a S0(3) matrix unlike the unit vector field in the NLSM for the Neel 
phase. The main feature of this new field theory is that there is no topological term 
and therefore no qualitative difference between integer and half- integer spin chains^^. 
There should be gap above the ground state(s) in either case. 

IV. RESULTS AND DISCUSSION 

We now present the phase diagram of the model and the various results obtained 
by us using DMRG. We also compare these with the existing analytical and numerical 
results. This is done in four subsections: A. The frustrated model with S — 0; B. The 
dimerized model with J2 = constant; C. The general J2 — S model; and, D. Coupled 
spin chains with 5 = 1. 

A. The frustrated spin chain ( J2 > 0, 5 = 0) 

Very few exact results are known for this model. The ground state is exactly 
solvable for J2 = 0.5. For cyclic boundary conditions, the ground state is doubly 
degenerate with the wave functions being given by the two possible Kekule structures, 
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namely, the two patterns of singlets which can be formed by neighboring sites. 

V^i = [1,2] [3,4] .... [2iV-l,2iV] 

(10) 

= [2iV,l] [2,3] .... [2Ar - 2, 2Ar - 1] , 
where [i, j] denotes the normalized singlet combination of the spins on sites i and j. 

The field theory studies by Haldane on the frustrated model revealed the existence 
of a transition from a gapless phase for J2 < J2c to a gapped phase for J2 > J2c ■ 
He estimated J^c to be 0.16. A direct computation of the gap obtained by exact 
diagonalization of small systems^ placed J2c near 0.30. However this turns out to be 
an unreliable method because of the essential singularity in Eq. (9) as J2 approaches 
J2C from above. The most reliable method of computing J2c is the one based on the 
crossing of the excited singlet and triplet states extrapolated to infinite system size 
through finite size scaling. Ref. 8 obtained a value of J2c = 0.2411 ± 0.0001 by this 
method. In the next subsection, we use this value of J2c to study the relation between 
the gap and 5. 

We have confirmed that the ground state is doubly degenerate for J2 > J2C ? 
as required by the Lieb-Schultz-Mattis theorem for gapped spin-half systems. Fig. 3 
summarizes our results for the extrapolated gap for < J2 < 1. The gap increases for 
J2 > J2c as a function of J2 until about J2 = 0.7 beyond which it decreases. In the 
limit of very large J2 , the gap should once again vanish since the system then consists 
of two decoupled nearest neighbor chains. We find that this behaviour sets in even at 
a relatively small value of J2 = 0.8. Fig. 4 shows the convergence of the gap with 
increasing chain length in the gapped ) as well as in the gapless {l/N) regions, 

and illustrates the validity of our extrapolation procedures. 

B. Dimerized model ( J2 = constant, 5 > 0) 

At J2c , we have calculated the gap A for various values of 6 on spin chains of 
upto 300 sites with the open boundary condition. The ground state is a singlet and 
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the lowest excitation is a triplet. The triplet is therefore obtained as the lowest energy 
state in the = 1 sector within the DMRG formalism. The minimum chain length at 
which the infinite chain behaviour can be expected to set in scales as the correlation 
length and hence depends inversely on A. For S — 0.007 (which is the smallest S we 
have studied), the infinite chain behaviour is expected for lengths greater than about 
200 sites. 

Fig. 5 shows that the log - log plot of the gap versus S is linear. The exponent 
for the gap is given by the slope of this straight line which is 0.667 ± 0.001. Previous 
numerical studies using exact diagonalization of spin chains with upto 22 sites had 
placed the value of this exponent between 0.9 and 1.0 I5,i6^ while those using finite size 
scaling placed it closer to 0.75^^. However these studies were all at J2 = 0, hence they 
were subject to errors due to the logarithmic corrections. Fig. 6 shows the log - log 
plot of the change in the ground state energy Eo{S) — Eo{0) versus 5. The slope of the 
straight line is 1.251 ± 0.001 which is somewhat difi^erent from the field theory result 
of 4/3. 

For a 'relativistic' massive theory which differs from a massless theory by a small 
perturbation (by small we mean that the 'velocity of light' is not changed), the mass 
gap A and the correlation length ^ should be inversely related to each other. In Fig. 
7, we plot the product ^A versus 6 for J2 = 0, and find that this does appear to be the 
case. 

C. The general J2 — 5 model 

The ground state of the model is known to be exactly solvable along the line 
2J2 + S = 1. The presence of the dimerization 5 lifts the degeneracy of the two ground 
states at J2 = 0.5. One of the Kekule states continues to be the ground state while the 
other is no longer an eigenstate of the Hamiltonian. Choosing the state 'ipi in Eq. (9) 
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as a trial wave function, we obtain a variational energy per site equal to 

Eo = -^(1 + 5). (10) 

Using the Rayleigh-Ritz variational principle, it is then easy to show that Eq saturates 
the lower bound to { H ) provided^ 2J2 + S — 1. Fig. 8 shows the variation of the gap 
along this line. 

The general model exhibits a gap in the excitation spectrum for all points (J2, S) 
not covered by the special cases discussed above. The ground state energies for various 
values of J2 and 5 are given in Table 1. 

To characterize the ground state further, we computed the equal time spin-spin 
correlation function at many points in the phase diagram. The system sizes chosen 
for these calculations were varied with J2 and 5 so as to reasonably approximate the 
infinite chain behaviour. We then computed the static structure factor S{q) defined 
as the Fourier transform of the correlation function. Some information about the type 
of long-range order (LRO), if any, can be obtained from the dependence of S{q) on q. 
The classical limit {S 00) predicts Neel order for J2 < 0.25 and a (coplanar) spiral 
order for J2 > 0.25. In other words, S(q) has a peak at qmax = tt for J2 < 0.25 and 
Qmax — cos~-^(— I/4J2) for J2 > 0.25. The periodicity of the ground state is I'njqmax ■ 

The quantum model has no LRO but exhibits a short-range order characterized 
by some qmax (and a finite correlation length if there is a gap in the spectrum). Earlier 
numerical studies^ of the S — 1/2 system revealed marked deviations from the large-^ 
(classical) results. They showed that qmax — ^ for J2 < 0.5 and qmax < tt for J2 > 0.5, 
with qmax approaching 7r/2 for large J2 . 

We have calculated S(q) for rings of sizes 100 and 150 at the points indicated in 
Fig. 9. The structure factor peaks at qmax = tt for all points along the J2 = line. 
On the (5 = line, S{q) peaks at qmax = tt for J2 < 0.5 and qmax < tt for J2 > 0.5 
with qmax approaching 7r/2 for very large J2 as shown in Fig. 10. In the full J2 — S 
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plane, S{q) shows a very interesting change of behaviour across the hne 2J2 + 5 — 1. 
For points to the left of the line, S{q) peaks at Qmax = while for points to the right 
of the line, S{q) peaks at Qmax < tt. However, for points close to the line on the right 
hand side, S{q) shows a very broad distribution near q — tt and we cannot identify a 
clear maximum in the plot (Fig. 11). 

The correlation length C is a minimum along the entire line 2J2 + 6 = 1, reflecting 
a highly disordered ground state. The main features in the S(q) plots seem to be in 
conformity with this behavior of ^. We can therefore call 2 J2 + 5 = 1 a disorder line. 
The disorder line separates the ground states with q-max = and qmax < as shown 
in Fig. 9. 

D. Coupled spin chains {S — 1) 

For 6 = 1.0 and J2 > 0, the model corresponds to two coupled spin chains as 
discussed earlier. The phase diagram of coupled spin-half chains was studied in Ref. 18. 
We have examined the dependence of the gap on the scaled strength of the interchain 
coupling Ji = 2/J2. This was calculated for a 2 x 100 ladder. In Fig. 12, we see 
a gap in the spectrum even for very small values of Ji as shown. Our results for the 
excitation gaps are in agreement with the studies of Barnes et al^^ which suggested 
the presence of a non-zero energy gap between the singlet ground state and the triplet 
excited states for all antiferromagnetic interchain couplings. 

The convergence of the ground state energy per site for this system is shown in 
Fig. 13. The fluctuation in the ground state energy with the total number of sites is 
because of the alternation between even and odd number of sites on a single chain. The 
amplitude of the fluctuation is large for weak interchain couplings and gets damped as 
the coupling between the chains increases. The energy per site when the intrachain and 
interchain exchange constants are equal is 0.5780. This value is intermediate between 
the ground state energy per site for the isotropic Heisenberg anitferromagnet on the 
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two-dimensional square lattice (—0.67) and the ground state energy per site on a chain 
(—0.4431). The advantage of our method is that it is still treated as a single chain 
problem and thus avoids the pit falls of DMRG in two dimensions. The accuracy of 
the method for the ladder is the same as that for the single chain. 

V. SUMMARY 

To conclude, we have studied the NN and NNN antiferromagnetic spin-half Heisen- 
berg chain with a dimerization, S, in the NN exchange. When S is zero, there is a critical 
value J2c = 0.2411 below which the system is gapless and above which the system is 
gapped. For the dimerized chain with the critical value of the NNN interaction, the 
energy gap above the ground state scales as S to the power 0.667 ± 0.001, while the 
change in the ground state energy Eo(6) — scales with the power 1.251 ± 0.001. In 
the J2 — S plane, we find a disorder line 2 J2 + 6 = 1 to the left of which the static 
structure factor S{q) peaks at Qmax — To the right of this line, Qmax gradually 
decreases from tt to 7r/2 for large J2 . For 6 = 1, the model corresponds to coupled 
chains as on a ladder. The system is gapped for all values of antiferromagnetic Ji . 
The energy per site is —0.57804 for Jj = 1.0 which lies in between the values for the 
one-dimensional and two-dimensional (square lattice) NN antiferromagnetic systems. 
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Table Captions 

1. Ground state energy per site Eg for various values of J2 and S. 

2. Ground state energy per site for coupled spin chains for various values of Ji . 
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J2 


5 




0.2411 


0.000 


-0.401866 


0.2411 


0.007 


-0.402681 


0.2411 


0.014 


-0.403766 


0.2411 


0.020 


-0.404844 


0.2411 


0.028 


-0.406430 


0.2411 


0.040 


-0.409051 


0.2411 


0.057 


-0.413132 


0.2411 


0.080 


-0.419154 


0.2411 


0.160 


-0.442862 


0.2411 


0.230 


-0.465728 


0.2411 


0.320 


-0.496844 


0.10 


0.0 


-0.42517 


0.20 


0.0 


-0.40885 


0.25 


0.0 


-0.40045 


0.30 


0.0 


-0.39284 


0.40 


0.0 


-0.38028 


0.50 


0.0 


-0.37500 


0.60 


0.0 


-0.38079 


0.70 


0.0 


-0.39711 


0.80 


0.0 


-0.42138 


1.00 


0.0 


-0.48565 


0.45 


0.07 


-0.40130 


0.48 


0.10 


-0.41281 


0.55 


0.10 


-0.41610 


0.15 


0.20 


-0.46329 


0.25 


0.35 


-0.50727 


0.40 


0.50 


-0.56611 


0.15 


0.60 


-0.60033 


0.20 


0.80 


-0.67613 


0.30 


0.80 


-0.67966 


0.48 


0.80 


-0.69256 



Table 1 
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Ji 


Eo 


0.020 


-0.44320 


0.040 


-0.44337 


0.066 


-0.44434 


0.133 


-0.44320 


0.200 


-0.44741 


1.000 


-0.57804 



Table 2 
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Figure Captions 

1. Heuristic phase diagram for the spin- half chain. The sohd hne along 5 = from 
J2 = to J2c is gapless; the rest of the diagram is gapped. The line 2 J2 + 5 = 1 
separates the Neel phase from the spiral phase. 

2. Schematic picture of the antiferromagnetic exchanges in the chain. 

3. Dependence of the gap A on J2 for 5 = 0. 

4. Convergence of the singlet-triplet gap A with 1/N for J2 = 0.5 (gapped system) 
and J2 = 0.2 (gapless) at 5 = 0. 

5. Log-log plot of the gap A versus 5 for J2 = 0.2411. 

6. Log-log plot of the change in the ground state energy Eo{5) — Eo(fS) versus 5 for 
J2 = 0.2411. 

7. Plot of the product of the correlation length ^ and the gap A versus 5 for J2 = 0. 

8. Gap A versus 5 along the line 2J2 + 5 = 1. 

9. Behaviour of the structure factor S{q) in the J2 — 5 plane. 

10. Plot of Qmax (i^ degrees) versus J2 for 5 = 0. 

11. S{q) versus q (in degrees) for various values of J2 and 5. 

12. Gap versus the interchain coupling Ji for coupled chains. 

13. Convergence of the ground state energy per spin Eq with system size for coupled 
chains with Jj = 0.04. 
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